Valley interference effects on a donor electron close to a Si/Si02 interface. 
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We analyze the effects of valley interference on the quantum control and manipulation of an 
electron bound to a donor close to a Si/Si02 interface as a function of the valley-orbit coupling at 
the interface. We find that, for finite valley-orbit coupling, the tunneling times involved in shuttling 
the electron between the donor and the interface oscillate with the interface/donor distance in much 
the same way as the exchange coupling oscillates with the interdonor distance. These oscillations 
disappear when the ground state at the interface is degenerate (corresponding to zero valley-orbit 
coupling). 
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I. INTRODUCTION 

Miniaturization of the traditional Si transistors is, af- 
ter decades, still following the time-defying Moore's law, 
with a 0.7 reduction in size every two years. Soon the ef- 
fects of dopants disorder in the transistors threshold volt- 
age will be importanti and, in a few years, their behavior 
will be dominated by quantum effects. These technical 
improvements can be used in the advancement of one of 
the promising candidates for the physical implementa- 
tion of a quantum computer: doped Si.-!^!^ Si provides 
with a scalable framework and very long spin coherence 
times, very important qualities for the actual realization 
of a quantum computer. 

The doped Si proposal for quantum computing (QC) 
relies on the control and manipulation over donor elec- 
trons which shuttle between a donor site and an interface 
with a barrier (Si02) which is a few Bohr radii away. 
This process is driven by an external electric field per- 
pendicular to the interface. For those relatively large 
distances compared to the lattice spacing, the problem 
may be treated within the effective mass approximation. 
Based on the single-valley effective mass approach, we 
have recently shown theoreticalljs^i^ that a donor-bound 
electron can be manipulated between the donor and the 
interface by suitably tuning the external electric field. 
The characteristic field for which donor ionization takes 
place is found to be a smoothly decreasing (inversely pro- 
portional) function of the distance between the donor and 
the interface.^ This result is in good agreement with ear- 
lier estimates based on tight-binding calculations where 
the peculiar conduction band of bulk Si with six degen- 
erate valleys in the (100) directions was included,^ in- 
dicating that valley interference effects may not play as 
strong a role in defining this characteristic field as it does, 
for example, in the behavior of the inter-donor exchange 
coupling.® However, other quantities of interest for the 
quantum control of donor electrons close to an interface 
may be affected by the multivalley structure of the Si 



conduction band, a question that deserves closer investi- 
gation. 

In this paper we generalize the effective mass treatment 
for electrons which shuttle between a donor site and an 
interface by considering a two-valley structure. This is 
appropriate for the interface bound states, since the six- 
fold conduction band degeneracy in bulk Si is partially 
lifted at the interface, where the two valleys perpendicu- 
lar to it become lower than the parallel valleys. Valley- 
orbit coupling at the interface leads to a non-degenerate 
ground state, where these two valleys contribute with 
equal weights. We confirm that the behavior of the char- 
acteristic field as a function of the distance of the donor 
to the interface is qualitatively similar to that obtained 
in the single- valley approximation. The quantitative esti- 
mates are also dependent on the magnitude of the valley- 
orbit coupling at the interface, a parameter which is quite 
sensitive to the interface type and quality!^ We have also 
studied the amplitude of the gap between the two low- 
est energy states at the characteristic field, a quantity 
related to the electron tunneling time. In this case, the 
single- valley estimates remain qualitatively reliable if the 
valley-orbit coupling at the interface may be neglected. 
However, in the presence of interface coupling, valley in- 
terference effects cause the gap to become an oscillatory 
function (at the lattice parameter scale) of the donor- 
interface distance, with quite small values of the gap at- 
tained at particular distances, leading to arbitrarily long 
tunneling times for donors at these positions. 

The outline of this paper is as follows. In Sec. |TT] we 
introduce our model and the formalism, while results are 
detailed in Sec. lIIII In Sec. lIVl we discuss the results and 
conclude. 



II. MODEL 

We consider a P donor in Si, a distance d from a (001) 
oxide interface, under a uniform electric field perpendic- 
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FIG. 1; (Color online) Double well potential in the z-direction 
(Coulombic donor potential plus triangular interface/electric 
field potential) and the lowest levels valley-orbit splittings in 
the case ^ IK^I for F ~ Fc- Note that the splittings 

depend both on the valley-orbit couplings and on the pinning 
phases zj and zd- 



orbit couplingjii Fj is the interface envelope function^ 



Ft 



f{z) X g{p) 



-a(2+(i)/2 



/3 



(4) 



12 ; ■ V"" 

with OL and /3 variational parameters. In the case of an 
isolated well, the pinning site r/ should be near the in- 
terface, where the envelope electronic density is maxi- 
mum and the potential is most attractive. Here we take 
r/ — —d (namely, at the interface); as shown below, the 
exact location of rj is not relevant for our results, as long 
as it is near the interface, far away from the donor site. 

This problem has been previously addressed in the sin- 
gle valley approach. In order to investigate multivalley 
effects in the simplest approximation, we assume that the 
donor ground state only involves the z and ~z valleys: 



ular to it (see Fig. [T|). The donor is at r = and the 
interface is at r = (0,0, —d). The ground state of the 
electron at the donor is writteniS 
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where Fj^ are envelope functions. In the presence of an 
interface (which is considered to produce an infinite bar- 
rier at 2; = — d), we take^ 



Fl, = N{z + d)e-VpV"^+^Vfc^ 



(2) 



with a and b variational parameters, and N a normal- 
ization constant. Here (/)p(r,r£)) = ^Bloch^^*''^ ^ 
M^(r)e*''^''^''"'"°^ are the Si band edge Bloch states {fi = 
1,2, . ..6), |k^| =ko = 27r0.85/as» with asi = 5.4 A the 
lattice parameter. The reference site r^j is irrelevant for 
the individual Bloch states, for which it only contributes 
with a phase. However, for the superposition state such 
as given in Eq. it represents a common pinning site 
of the Bloch states which, for the ground state of an iso- 
lated donor in Si, is naturally chosen to be at the position 
of the donor r^j = OJ^ 

At a (001) interface, the conduction band edge six-fold 
degeneracy is broke n^^'^^ *^ and the z and —z valleys 
are ~ 40 meV lower in energy than the other four val- 
leys. The lowest energy states at the interface are the 
symmetric and antisymmetric combinations of these two 
valleys 



*f = ^^/(r)[0.(r,r,)±0_,(r,r,)] 
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z and —z are the relevant components which couple more 
strongly to the interface lowest states, and the lowest 
ones in Si under tensile strain [namely, grown on relaxed 
Sii_a;Gea; (001)]. ~ The envelopes Fu and Fj are the same 
as obtained in the single- valley calculation where the val- 
ues of the variational parameters a, 6, a, and /3 were 
determined. The remaining parameters, r^i = (0,0, z/j) 
and rj = (0, 0, z/), are optimized variationally here. 

The Hamiltonian of a donor electron close to a Si/Si02 
interface is, in rescaled atomic units a* = h^esi/Ta^e^ — 
3.157 nm and Ry* = m^e^/2h^el^ = 19.98 meV, 



H 



2 

KcFz — 



■ 7 



dz'^ 
2Q 



Vp^TTITW 2(z + d) 



(6) 



where 7 = m^/m<^\, Q = (esiOs - esi)/(esi02 + esi), 
esi = 11.4 and esiOs = 3.8, k = 3.89 x 10^'^e|, {m/m^f 
cm/kV, and the electric field F is given in kV/cm. The 
last term describes valley-orbit effects, namely the cou- 
pling between different valleys due to the singular nature 
of both the donor [D) and the interface (/) potentials. 
These couplings are quantified by the parameters and 
1^0, as described below. 

In our two-valley formalism, the problem is restricted 
to the basis set of the lowest uncoupled donor and inter- 



im) face states, namely {Fu"^ 
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whose degeneracy is lifted by the abrupt interface valley- Hamiltonian matrix 



Bloch ■= 



Bloch 



^/^Bloche 



-ikozo 
ika{zi+d) 



}, leading to the 



3 



H 



/ Ed 

Vj^ COsi2koZD) 

V 



V;^ cos(2fco2D) Hioe'''"'^'^-"'^ 
Ed 
Ei 





Vicos[2ko{zi + d)] 
Ei 



(7) 



where Ed = (Fc|i?|i^z3), Ej = and = 

are the same as the single valley matrix ele- 
ments. The two additional parameters, the pinning posi- 
tions zd and zj, are calculated variationally. The valley- 
orbit (VO) coupling parameter at the donor is known 
from the isolated donor spectrum (splitting within the 
IS manifold): For a P donor it is = -1.5 meV^ii^ 
We use this value throughout. However, the VO coupling 
at the interface depends on the sample, and different val- 
ues have been reported in the literature!^ We take V^^ 
as a parameter varying from to —10 meV. The nega- 
tive sign is chosen in analogy with the VO coupling at 
the donor, but our main conclusions only depend on the 



absolute value of V^q. For Vyo < 0, the ground state 
is a symmetric combination of the z and —z valleys. 
Full valley-orbit coupling (maximum splitting between 
the symmetric and antisymmetric combinations of the z 
and —z valleys, given in Eq. ([3]), with the symmetric one 
as the ground state) corresponds to zd = at the donor 
well and zj ~ — d at the interface well. When V^^ — 0, the 
two combinations are degenerate at the interface. This 
degeneracy is broken when Vj^ cos[2ko{zj + d)] ^ 0. Note 
that the donor states, defined in Eq. [5l would also be de- 
generate if ZD is such that cos(2A:o2:d) = 0. 

We define the overlap matrix 
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where S — {Fd\Fi) is the overlap between the interface 
and donor envelopes, which is a decreasing function of d. 
The four eigenvalues Ei and eigenstates Xi {i = 0, 1, 2, 3) 
satisfy Hxi = EiSxi. 

According to our model, meaningful calculations refer 
to electric field intensities for which the envelope func- 
tion of the donor well ground state may be described by 
a single IS- like function. As shown in calculations of the 
Stark effect of P in Sij^ higher donor excited states con- 
tribute to the ground state under fields above 24 kV/cm. 
We thus restrict our studies to intermediate to long dis- 
tances, namely d > 5.5a*, so that the characteristic fields 
(see Fig. [1]) remain below ~ 24 kV/cm. 



III. RESULTS 

As discussed in previous studies,® the characteristic 
field Fc above which the electronic ground state is at 
the interface, as well as the time required to 'shuttle' the 
electron from the donor to the interface, are determined 
by a level anticrossing process involving the electronic 
ground state. The tunneling time r has been estimated 
from the gap g at anticrossing, t = fi/g. We focus here 
in the study of Fc and g as a function of d and of y^Q, 



the interface VO coupling parameter. 

In Fig. [T] we illustrate the double well potential (donor 
plus interface) and the uncoupled two- valley levels (i.e., 
assuming Hjd = 0). In the single valley approxima- 
tion,? the problem is two-dimensional, only two levels 
(one per well) are relevant and one anticrossing occurs at 
Fc- In the present case, for V.^^ ^ 0, there are two non- 
degenerate levels in each well and, therefore, there are 
four level anticrossings as illustrated for V^^, = ~1 meV 
and d = 5.94a* in Fig. [2] (a). The relevant anticrossing to 
determine Fc and the gap g is the one that involves the 
lowest eigenstate. When the interface states are doubly 
degenerate, namely V^^, — ^^ there are only two level an- 
ticrossings involving three levels each, see Fig.[2](b). One 
of the levels at the anticrossing remains uncoupled to the 
other two that present the typical level repulsion with an 
associated gap. For < \V^g\ < 0.02 meV, the results are 
sensitive to d though the predominant behavior is similar 
to the degenerate case shown in Fig. [2] (b). 

In Fig. [3] we compare the values of the characteristic 
field Fc required to move the donor electron ground state 
from the donor to the interface within the single valley 
approximation and the 2-valley formalism for V^^^ = 
and V^c, = ~1 meV. The behavior of Fc obtained in the 
2-valley formalism is qualitatively similar to the single 
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FIG. 2: (Color online) (a) Here we show the eigenvalues for 
d = 5.94a* and V^^ = ~1 meV. All the levels are coupled and 
there are anticrossings between the levels in pairs. The insets 
are blow-ups of each of the four level anticrossings. The rel- 
evant anticrossing is the one involving the ground state. We 
get qualitatively similar results for \Vjo\ ^ 0.02 meV and for 
all distances d not at a gap minimum (see Fig. |4] (a)), (b) 
Levels dependence on the electric field when the interface VO 
coupling is exactly zero Vjo = 0. We show the eigenvalues for 
d = 5.99a* . The plot is qualitatively the same for all distances 
d. The insets are blow ups of each of the level anticrossings 
which shows three lines (an intermediate level remains uncou- 
pled to other two, which display the anticrossing structure) 
due to the double degeneracy of the interface states. 



valley approximation; only a small shift is obtained due 
to the different values of the ground state levels at each 
of the wells. The shift with respect to the single valley 
result is dependent on the relative values of \V^J ^^"^ 
\V^\, being zero if \VjJ — \Vj^\. For a particular value 
of d, Fc decreases as \VjJ increases. 

The smooth behavior of Fc in the 2-valley formalism 
is not shared by the other quantity of interest: the gap 
at anticrossing (related to 'shuttling' tunneling times). 
This is illustrated in Fig. 2] (a) where oscillations for the 
case T^o = ~1 meV are patent. The maxima of the gap 
oscillations correspond to tunneling times of the order of 
the single valley results^ [see Fig. 0] (b)] while the gap 
minima would lead to much longer tunneling times. The 
periodicity of the oscillations is given by abs[cos(fco(i)] 
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FIG. 3: (Color online) Fc versus d for the single valley ap- 
proximation, V^o = ~1 meV, and Vj^ ~ (bottom to top 
curves). No oscillatory behavior is obtained as expected from 
the tight-binding results in Ref. {Tj. The shift between the dif- 
ferent curves is due to the relative values of IV^^"^! (fixed to the 
bulk value Vj^ = —1.5meV) and \Vja\- 



which is not commensurate with the lattice. Note that 
d is not a continuous variable in a real system but its 
value changes in monolayer steps of 1.35A. These values 
of d correspond to the stars that sprinkle the continuous 
oscillating dashed line in Fig. 2] (a). Therefore, tunneling 
times are in practice undetermined (and could be much 
longer than predicted in the single valley approximation) 
unless the position of the donor is known with atomic pre- 
cision. The gap amplitude (defined by a phenomenolog- 
ical fitting discussed in the figure caption) also depends 
on the value of V^^ as shown in Fig. |4] (b). Increasing 
\Vj^\ leads to smaller gaps at anticrossing and therefore 
longer tunneling times (the difference between maximum 
and minimum times at a particular d is ^ 20%). In con- 
trast, for — 0, the gap decreases smoothly (with no 
oscillations) as d increases. 

The qualitatively distinct results obtained for V^^ = 
and V^^ — ~1 meV illustrated in Fig. |4] (a) emerges from 
the distinct level structures for these two cases, as shown 
in Fig. [21 For Vj^ = 0, the purely interface states are 
degenerate (high and low field limits in this figure). Un- 
der particular fields, this degeneracy is lifted because a 
specific superposition couples to a donor state, leading 
to an anti-crossing between these levels, while a second 
interface state remains uncoupled and unaffected. The 
relevant anticrossing in our studies does not involve the 
uncoupled state, which is never the lowest (or highest) 
energy state within the anticrossing range of fields. For 
0, our numerical investigations show that for very 
small values of this parameter (e.g. \VjJ ~ 0.01 meV) 
and for some values of d, an interface level remains uncou- 
pled and the gap behavior is very similar to the V^^ — 
case. The oscillations in the gap appear consistently for 
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FIG. 4: (Color online) (a) Gap versus d for V^^ = ~1 meV 
(stars and oscillating thin dashed line), V^o = (dashed line) 
and single valley approximation (solid line). The thin con- 
tinuous dashed line is a fitting of the form abs[cos(fc()d)] x 
[Aexp{-d/B)] with A = 5.03248i?y*, and B = 0.82896a*. 
The exponential is a fitting to the curve formed by the max- 
ima of the oscillating function. The stars correspond to mono- 
layer steps of I.35A. The curve for Vjo = can be also fitted 
to ^exp(-d/B) with A = 4.98483i??/* and the same B. (b) 
Value of the fitting parameter A (or 'gap' at d = 0) as a 
function of V^o (in meV). The fitting parameter B is always 
B = 0.82896a*. Changing the sign of V^^ leads to a shift of 
half a period in the gap oscillations while A{V^o) = A{—Vyo)- 



all d and \Vj^\ > 0.02 meV. 

In Fig. [5] we show the values of the phase pinning sites 
zj and Z]j obtained variationally. The pinning sites de- 
pend on the applied electric field. Note that the splitting 
between the eigenstates in each of the wells depends on 
the cosines cos{2koZ£)) and cos[2fco(z/ -I- d)] (see Fig. [T|). 
Therefore, in a generic case as the one depicted in Fig. [5] 
(a), the splitting between the two eigenstates (the sym- 
metric and antisymmetric combinations of the z and —z 
valleys) depends on F. Importantly, the splitting is max- 
imum [cos{2koZD) = ±1 or cos[2A:o(z/ + d)] = ±1] in the 
well where the electron ground state is (namely, at the 



FIG. 5: (Color online) Optimized variational parameters zi 
and ZD, related to the pimnng phases, for Vjo = —1 meV. 
This figure illustrates how the oscillations in the gap emerge 
from the phases in the cosines in Eq. (O . (a) Cosines squared 
of the pinning phases at d = 6.2a*. The two curves cross 
at Fc- The pinning shifts from the donor at small fields 
(namely, cos^(2fco2D) ~ 1) to the interface at larger fields 
(namely, cos^ [2^0(2:/ + d)] — 1). The value of the cosine 
squared in the unoccupied well depends on d and oscillates 
between and 1. For instance, in the case depicted in (a) 
when cos^(2fco2D) ~ 1, cos^(2A:o(2:/ -I- d)) goes down to ~ 0.2 
(point signaled with an arrow at the lower left corner). In (b) 
we plot the value of cos'^(2fco(zj-l-d)) at small fields {F <^ Fc) 
for different d's (stars) and we see that it follows abs[cos(fcod)] 
(continuous line), cos^ (2^020) shows the same behavior when 
F ^ Fc- These oscillations are the same shown by the gap, 
illustrated in Fig. |3] (a). Therefore, for a d corresponding 
to a gap minimum, the cosines squared go from to 1 (and 
viceversa), while for a d corresponding to a gap maximum, 
both cosines remain very close to 1 for all electric fields. This 
means that for a d corresponding to a gap minimum, the pin- 
ning phase shifts from donor to interface with increasing elec- 
tric field while for the cases corresponding to a gap maximum 
each state remains pinned at its origin (donor and interface 
respectively) for all fields. 
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donor for F < Fc, and at the interface for F > Fc). The 
sphtting in the weU where the excited state is (the inter- 
face aX F < Fc and the donor a.t F > Fc) oscillates with 
d between zero and maximum with the same periodicity 
as the gap, as shown in Fig. [5] (b). This implies that for 
the values of d close to the minimum gap (cos(fcoc?) ^ 0), 
the eigenstates at the donor (interface) are almost degen- 
erate for F > Fc {F < Fc). In these cases, the splitting 
switches from zero to its maximum possible value, and 
viceversa, at Fc, with anticrossings involving the 4 lev- 
els and producing a minimum gap. When d corresponds 
exactly to degenerate states [cos{kod) — 0], the gap is 
larger than the ones depicted in Fig. [H However, due to 
the incommensurability of cos(kod) with respect to the 
lattice, the chance that d satisfies cos{kQd) = is very 
small. 

If zj and zo are fixed such that the related cosines 
are finite at all fields, the oscillations in the gap are pre- 
served. The oscillations disappear when at least one of 
the wells has a degenerate ground state at all fields, and 
the problem reduces to the single valley approach if both 
VO couplings are zero. 

IV. DISCUSSIONS AND SUMMARY 

It has been known for a few years that the exchange be- 
tween electrons bound to donors separated by a distance 
R oscillates with R due to the interference between the 
valleys in the conduction band of SijSiAZ These oscilla- 
tions present consequences for Si quantum computing if 
the 2-qubit operations rely on exchange coupling. Here 
we find that valley interference is also relevant for the 
quantum control and manipulation of single electrons in 
a donor/interface system. In particular, the tunneling 
times (or gaps) for the donor-interface shuttling have 
a periodic (non-commensurate with the lattice) depen- 
dence on the distance d between donor and interface. 
These oscillations arise only when Vjc, is non-zero {V^ 
is fixed to the bulk value of —1.5 meV). For Vjc, = the 
qualitative behavior of the gap with d reduces to that 
of the single valley approximation, namely, it shows no 
oscillations. 

The oscillations in the gap versus d can lead to very 
long tunneling times for donors at distances such that 
cos(fcod) ~ 0. In the current experiments, donors are 
ion implanted and form a doped layer with a distribu- 
tion of distances with respect to the interface!^ In this 
case, our results would imply that some of those donors 
would never ionize in the relevant time scales, producing 
an effective reduction in the planar density of donors, 
a desirable feature in terms of allowing experiments to 
approach the single donor behaviori^ 

On the other hand, arbitrarily long tunneling times for 
donors at certain distances represent a serious overhead 



for the final objective of implementing a quantum com- 
puter using doped Si, as we need to be able to manipulate 
the electrons many times within the spin coherence time 
(> 1 ms in bulki^iSSii^) . The disappearance of the oscil- 
lations when Vjci ~ is of no help as a finite valley-orbit 
coupling is required for a well defined qubit, in particular, 
the valley-orbit splitting has to be larger than the spin 
splitting, which would determine the two states of a spin 
qubit. Fortunately, the statistical weight of the positions 
corresponding to vanishing gaps is relatively small (see 
the stars in Fig. 2] (a)) so most donors will be at positions 
at which the tunneling times are of the same order as the 
ones estimated within the single-valley approachi^ 

Experimentally, different values of Vj^ have been re- 
ported depending on the quality of the interface and 
growing method.^ This coupling can be a complex num- 
ber. We have chosen V^^ to be negative and real and 
explored a wide range of values > V^^, — ~10 meV 
consistent with the ones reportedi^ We have also checked 
that, for Vjg > 0, the gap versus d amplitude is coin- 
cident with the case Vj^ < though its oscillations are 
shifted by half a period. We expect that a complex Vjc, 
would give different shifts in the periodic function with- 
out affecting the gap amplitude and the decaying behav- 
ior (given by the fitting parameters A and B in Fig. [J]). 

In conclusion, we have analyzed the effect of the multi- 
valley structure of the Si conduction band on the manip- 
ulation of electrons bound to donors close to a Si/Si02 
interface. Our results are based on a two-valley effec- 
tive mass formalism, introduced here, and our main con- 
clusions are expected to remain valid if the full multi- 
valley degeneracy at the donor (with six valleys) were 
included. The characteristic field needed to shuttle an 
electron between the donor and the interface well is not 
qualitatively affected by valley-interference. The behav- 
ior of this quantity is in agreement with our previous 
single-valley treatmentjS However, the tunneling times 
involved in the process oscillate with the distance be- 
tween the donor and the interface at the atomic scale. 
The single- valley results typically give a lower bound for 
the tunneling times (upper bound for the anticrossing 
gaps). This behavior must be taken into consideration 
in the practical implementation of a doped Si quantum 
computer, as direct and full control over the operations 
discussed here may require positioning of donors within 
atomic precision. 
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